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Forward 


Data analysis by AI is a common thing 
nowadays. In the field of data analy- 
sis, many problems are being solved 
by machines by applying this tech- 
nique.Especially Deep Learning tech- 
nique has gained prominence in the 
time of now. Medical Science, Agri- 
culture, and Bio-Informatics are the 
free movement of AI and Deep learn- 
ing in all fields. In the field of combat 
modeling is no exception. Driven by 
this, UAVs, UGVs, and drones have 
now taken a special place in the era 
of Autonomous forces.To make these 
autonomous devices more sophisticated, 
they need to train more with histor- 
ical battle data. As Human Soldier 


ii 


needs to be trained by Historical Bat- 
tle Studies similarly Autonomous Robotic 
Forces need to be trained by histor- 
ical data.This book is written very 
nicely to explain how battle data can 
be used for train the autonomous forces 
with the Sun Tzu’s Art of War phi- 
losophy : 

Happy reading! 


Prof. Guang Yang 


ICL,London,UK 
Date:3"“Feb,2024 
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Preface 


Recent times RL has gained signif- 
icant popularity in simulating com- 
petitive games. The ability to re- 
inforcement learn patterns through 
deep artificial neural network is an 
essential components of an intelligent 
force.: 


1. Deep learning is therefore an 
ideal choice of so called "Intelli- 
gent military robots" which in- 
volve processing and fusion of 
data from different sources and 
historical databases of differen- 
tial equation based combat mod- 
elingCombat modeling. 


2. Advance level of Artificial In- 
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telligence (AI) based Model Fit- 
ting 


3. Introduction of simple theorems, 
Notations, and logic of combat 
modeling. 


4. The glossary and Acronyms of 
most useful terms of combat mod- 
eling. 


This book is a beginner’s introduc- 
tion to Artificial Intelligence based 
model fitting for generating simple 
equations which can also be used in 
combat analysis, defence planning or 
even on in any battle analysis. The 
book does not attempt to compile all 
modeling techniques and algorithm 
available in the open literature it tries 


to present those models which suits 
best with combat analysis. The goal 
is to keep all important historical pat- 
tern setting battles in a common place 
which can be shared as graphical and 
network presentations in the defense 
analyst. 

Happy reading! 


Sumanta Kumar Das 


New Delhi 
Date:3"“Feb,2024 
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Summary 


Can Autonomous Forces be trained 
with Artificially Intelligence (AT) al- 
gorithms like Deep learning from his- 
torical data and further be improved 
by the rules set by Sun Tzu’s Art 
of War for developing a state-of-the- 
art system of planning and analysis 
of Military Strategy. The applica- 
tion of Machine learning (ML) algo- 
rithms are general practice for ex- 
tracting meaningful models from the 
data.ML could be unsupervised, su- 
pervised, or reinforcement learning. 
This book explores deep learning neu- 
ral networks for extracting model pa- 
rameters for generating autonomous 
forces, Although our main focus is on 
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the world's largest tank battle that 
is Battle of Kursk, we are also ex- 
ploring other historical battle data 
of the world for validating the sys- 
tem. The battle of Kursk between 
the Soviets and Germans is known 
to be the biggest tank battle in his- 
tory. This book explores the two- 
dimensional-simplex tank and artillery 
data from the Kursk database for an- 
alyzing a class of discrete time homo- 
geneous and heterogeneous Lanchester 
models. Under homogeneous form, 
the Soviet’s (or German’s) tank ca- 
sualty is attributed to only the Ger- 
man’s (or Soviet’s) tank engagement. 
For heterogeneous form, the tank ca- 
sualty is attributed to both tank and 
artillery engagements. For validat- 
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ing the models, different. goodness- 
of-fit statistics are used for compari- 
son. 
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Abstract 


Can autonomous forces 
be trained with Artificial In- 
telligence (AI) algorithms like 
deep learning from histori- 
cal data and further be im- 
proved by the rules set by 
Sun Tzu’s Art of War for 
developing a state-of-the-art 
system of planning and anal- 
ysis of military strategies? 
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The application of Machine 
Learning (ML) algorithms 
are general practice for ex- 
tracting meaningful models 
from the data.ML could be 
unsupervised, supervised, or 
reinforcement learning. This 
book explores deep learn- 
ing neural networks for ex- 
tracting model parameters 
for generating autonomous 
forces, Although our main 
focus is on the world’s largest 
tank battle that is Battle of 
Kursk, we are also explor- 
ing other historical battle 
data of the world for vali- 
dating the system. The bat- 
tle of Kursk between the So- 
viets and Germans is known 
to be the biggest tank bat- 


tle in history. This book 
explores the two-dimensional- 
simplex tank and artillery 
data from the Kursk database 
for analyzing a class of dis- 
crete time-homogeneous and 
heterogeneous Lanchester mod- 
els. Under homogeneous form, 
the Soviet’s (or German’s) 
tank casualty is attributed 
to only the German’s (or So- 
viet’s) tank engagement. For 
heterogeneous form, the tank 
casualty is attributed to both 
tank and artillery engage- 
ments. For validating the 
models, different goodness- 
of-fit statistics are used for 
comparison. 


Keywords— Battle of Kursk, Ho- 


mogeneous , Heterogeneous Lanch- 
ester Model, Goodness-of-fit, Sum- 
of-square residuals, Chi-square, Kol- 
mogorov Smirnov. 


1 Introduction 


During the past decades, many dif- 
ferential equation based models [20, 
26, 32, 43, 36] have gained significant 
importance for representing combat 
dynamics. These equations are widely 
used for modelling in warfare and rep- 
resenting the decrease in force levels 
over time commonly referred to as 
attrition process.Lanchester in 1914 
proposed a set of differential equa- 
tions, which quantify the importance 
of force concentration on the battle- 
field. Many authors have subsequently 
modified his original work to repre- 
sent combat dynamics in modern war- 
fare. In the recent time Al-enabled 
warfare involving "Killer Robots" [2] 


has forced the designer to change the 
visualization of combat model devel- 
opment from traditional statistical meth- 
ods to more practical adaptive and 
heuristic approaches. In [10] the au- 
thors have discussed how artificial in- 
telligence (AI) could be used in political- 
military modeling, simulation, and wargam- 
ing of conflicts with nations having 
weapons of mass destruction and other 
high-end capabilities involving space, 
cyberspace, and long-range precision 
weapons. In reference [6], we found 
that a fast large-scale theater model 
combining ground to ground battle 
attrition with air-to-ground strikes has 
been developed using such models. The 
features of the Lanchester’s Equation 
that makes it suitable for analysis in- 


cludes: 


Applicability: Lanchester mod- 
els are widely used for histori- 
cal battle analysis [14]. Other 
than analysing human warfare 
Lanchester model have also been 
used for analysis of fights among 
social animals, market analysis 


[1]. 


Force Aggregation: Lanchester 
models are found to be suitable 
for developing aggregated com- 
bat modelling using High Res- 
olution Simulation Model. In 
reality, actual historical com- 
bat data is not easily available 
and common practice is to de- 
velop High Resolution simula- 


tion data with detailed design. 
Various literature [5], [9] [35] 
have demonstrated that estimat- 
ing attrition rates from high- 
resolution simulation and using 
Lanchester model for linking the 
various resolution of different 
simulation model. 


e Flexibility: Lanchester mod- 
els are flexible for both homo- 
geneous as well as heterogeneous 
situations.Lanchester models are 
used for theoretically consistent 
force aggregation and dis-aggregation 
in two dimensions [38, 21]. 


Regardless of credits of prior discov- 
ery, Lanchester's equations are used 
worldwide for calculating attrition rates. 


We propose a general form of the het- 
erogeneous Lanchester's model as: 


F 


i=l 


F 
Y = AVE XE Vi iiu 


i=1 
(2) 

where X;, denotes the strength of the 
i th type of Red forces at time t and 
Y;, denotes the strength of the ith 
type of Blue forces at time t. Xi 
and Y,, are red and blue forces killed 
at time t. 

a; represents attrition rate of ith 
type of Blue forces and 

b; represents attrition rate of ith 
type of Red forces; 
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Vi —1,2,..F 

where 

F denotes the total number of forces. 

pi the exponent parameter of the 
attacking force, 

qi is the exponent parameter of 
the defending force. 

Equations (1) and (2) involve un- 
known parameters a;,b; ,p; and q;. A 
lot of work has been done to estimate 
these parameters using statistical es- 
timation methods like Least Square, 
Maximum Likelihood, Bayes, Method 
of Moments [9, 35] etc. These es- 
timates are most suitable for homo- 
geneous situation. In heterogeneous 
situations these estimation procedures 
fail. In this book we propose AI based 
generalized reduced gradient (GRG) 
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algorithm |18, 30| which is alterna- 
tively solved through Deep learning 
Neural Network. We are generally 
acquainted with two forms of these 
equations for homogeneous weapon 
engagement (when i. Lanchester lin- 
ear law in which p; = qi = 1 and 
force ratios remain equal if a XP = 
b;.Y;"". Lanchester's linear law is in- 
terpreted as a model from a series 
of one-on-one duel between homoge- 
neous forces and this law describes 
combat under "ancient conditions". 
The equation is also considered a good 
model for area fire weapons, such as 
artillery. Lanchester square law in 
which , that is, force ratios remain 
equal if applied to modern warfare 
in which both sides are able to aim 
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their fire or concentrate forces. 
On integrating equation (1) and 
(2) we obtain the state equation: 


XN Yn 
> VAINA 
w=1j=1 
XN Yn 
IAN 


i=1j= 


=1 (3) 


where X?, X9, Y, Y? , represent the 
initial values of Blue and Red forces 
respectively. This equation says that 
the relationship between the power 
of the losses in any fixed time period 
is equal to the inverse ratio of the 
attrition rate parameters. Equation 
(3) leads to the victory condition for 
Blue. Most forces have breakpoints 
at which they will cease fighting and 
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either withdraw or surrender if: 


Xn YN . 
2; 2, ai X; X; YiYj 


i=1 j=1 


UTR >1 (4) 
X O 
i=1 j=1 


Finally equations (4) may be solved 
in closed form as function of t. 


i 1 . . Qi a; b. 
XP = (raa em 


1 , ; Qi, t /a;.b; 
eat tee 
(5) 
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. 1 ! | dbi ab. 
Y? = (vg - xP eve) 


a ai 


shee — xb. yea, 
(6) 
There is another form of mixed 
combat model where attacker uses area 
fire (pj = q; = lie. linear form) against 
a defender using aimed fire ( p; = 1,q; = 0 
i.e. square form). This mixed form 
of combat model is known as ambush 
model proposed by Deitchman [11]. 
Helmbold [19]in 1965 studied the 
Iwo-Jima campaign between USA and 
Japan using one-sided homogeneous 
Lanchester model. Bracken in 1995 
studied Ardennes Campaign between 
Germany and USA. Clemens in 1997 
[7] and Lucas and Turkes in 2003 [29] 
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studied the Kursk campaign between 
Soviet and Germany. Willard [46] 
has tested the capability of the Lanch- 
ester model for analyzing the histor- 
ical battle data for the battles fought 
between the years of 1618-1905. Bracken[4| 
(1995) used the database of the Ar- 
dennes campaign of World War II for- 
mulating four different models which 
are the variations of the basic Lanch- 
ester equations. The models devel- 
oped in his study were homogeneous 
in nature in terms of tank, APC, ar- 
tillery etc. He concluded that Lanch- 
ester linear model best fits the Ar- 
dennes campaign data in terms of min- 
imizing the sum of squared residuals 
(SSR). This work validates the ap- 
plicability of the Lanchester model 
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for the historical Battle data. Fricker 
[16] revised the Bracken’s models of 
the Ardennes campaign of World War 
II. He extended Bracken’s model by 
applying linear regression on the log- 
arithmic transformed Lanchester equa- 
tions and included the data from the 
entire campaign and air sortie data 
as well. Lastly, he concluded that 
neither of the Lanchester linear or 
square laws fit the data. A new form 
of Lanchester equations emerges with 
a physical interpretation. 

Clemens [7] fits the homogeneous 
version of Lanchester equations to the 
Battle of Kursk. He used two differ- 
ent techniques (i) Linear regression 
on logarithmic transformed equations 
(ii) a non-linear fit to the original 
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equations using a numerical Newton- 
Raphson algorithm. 

Hartley and Helmbold [17] exam- 
ined the validity of Lanchester’s square 
law using the one-sided data from the 
Inchon-Seoul Campaign. They have 
not found good fit using constant co- 
efficient square law but better fit was 
found when the data was divided into 
a set of three separate battles. They 
concluded the Lanchester’s square law 
is not a proven attrition algorithm 
for warfare although they also com- 
mented that one-sided data is not suf- 
ficient to verify or validate Lanch- 
ester square law or any other attri- 
tion law. They have used linear re- 
gression, Akaike Info criterion and 
Bozdogan’s consistency AIC(BAIC). 
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Based on the regression analysis they 
have found the models with three re- 
gression parameters with intercept and 
without intercept was the best model 
with higher degree Coefficients of de- 
termination. 

NR Johnson and Mackey [22]anal- 
ysed the Battle of Britain using the 
Lanchester model. This was a bat- 
tle of an air combat between German 
and Britain. 

Wiper, Pettit and Young [47] ap- 
plied Bayesian computational tech- 
niques to fit the Ardennes Campaign 
data. They studied stochastic form 


there is role of any attacking and de- 
fending army on the number of casu- 
alties of the battle. They compared 
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of Lanchester model and enquired whether 


their results with the results of the 
Bracken and Fricker and results were 
found to be different. They concluded 
that logarithmic and linear-logarithmic 
forms fits more appropriately as com- 
pared to the linear form found by 
Bracken. They also concluded that 
the Bayesian approach is more ap- 
propriate to make inferences for bat- 
tles in progress as it uses the prior 
information from experts or previous 
battles. They have applied the Gibbs 
sampling approach along with Monte 
Carlo simulation for deriving the dis- 
tribution patterns of the parameters 
involved. 

Turkes [40] extended the previous 
work for the validation of Lanchester 
models with real data. He stated 
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that historical data for validation of 
attrition model is poor. Mostly, the 
data contained starting sizes and ca- 
sualties only for one side. He applied 
various derivatives of Lanchester equa- 
tions for fitting model on the Kursk 
Database. The results found in his 
study were different with earlier stud- 
ies on the Ardennes campaign. He 
found that wide variety of models fit 
the data as well. He has shown none 
of the basic Lanchester models fit the 
data, bringing into question their use 
in combat modelling. Lucas and Turkes 
[29] used a new approach to find the 
optimal parameters for fitting Lanch- 
ester models on the data of Battles 
of Kursk and Ardennes. They have 
gained an understanding of how well 
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various parameter combinations ex- 
plain the battles. They have found 
that variety of models fits the data. They 
concluded that none of the basic laws 
(i.e. square, linear and logarithmic) 
fit the data correctly and raises the 
question of utility of basic Lanchester 
model for combat modelling. They 
also suggested finding new ways to 
model the aggregated attrition pro- 
cess to provide a good-fitting Lanch- 
ester model. 

US Army’s Colonel Trevor N. Dupuy 
[12] discussed about different weapon 
system in his book about the Evolu- 
tion of Weapons and Warfare (1990) 
which had evolved from 2000 BCE 
onwards till the Cold War and their 
tactical impact on combat. Despite 
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its Western bias, the book is good for 
detailed description of the military 
hardware which modern Europe pro- 
duced. Eminent author K Roy [34] 
described a global history of warfare 
from slings to drones also includes 
discussion on insurgency, civil war, 
sleges, skirmishes, ambushes and raids. 
The main aim of this book is to 
fit Lanchester Model based on Kursk 
data.For that we require to estimate 
attrition rates and exponent param- 
eters. There are several approaches 
to estimate the parameters. We shall 
consider two common and rational 
procedures namely, Least Square Es- 
timation (LSE) and Maximum Like- 
lihood Estimation (MLE). These two 
estimation procedures will be applied 
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by AI techniques like Deep Learn- 
ing. The authors in [27] have demon- 
strated Deep Reinforcement Learn- 
ing (DRL) AI techniques for solving 
complex problem, our autonomous forces 
are configured with this technique. 
DRL is an iterative optimization pro- 
cess. Our autonomous force higher 
commanders are established on the 
logic of Deep Deterministic Policy Gra- 
dient (DD PG Agent) [27]. These 
DDPG Agents are represented by actor- 
critic relationship where Actor is (S, 0), 
for observations S with parameter 0 
and critic Q(¢,S,A) with parameter, 
observation and action are ¢,S,A re- 
spectively. The parameter 0 is peri- 
odically updated for optimal value of 
long term reward function that pro- 
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duces the target Actor $(5,0,) and 
target critic Q(¢;,5,A). In the present 
book the observations are historical 
battle data, actions are mathemati- 
cal functions formulated from "Art of 
the War"[41],the parameters 0 esti- 
mated are attrition rate coefficients, 
exponents parameters which are be- 
ing estimated through Deep Learn- 
ing Neural Network, the critic pa- 
rameters $, are estimated from the 
same data set with strategic inputs 
and rules from the "Art of the War". 
'The long term reward function is de- 
signed from the GOF measures as de- 
fined in subsequent sections. The ob- 
servations sets are Soviet’s and Ger- 
man's Tank, APC, Arty Gun , In- 
fantry. That is 1 * 8, these datasets 
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are discrete in nature. The actions 
are casualties in these components 
which is recorded as discrete elements. 
To create the discrete and continuous 
observation and action spaces for the 
RL, RINumericSpace and RlFinite- 
SetSpace are used. 

In the next section we have dis- 
cussed in detail the mathematical for- 
mulations of homogeneous and het- 
erogeneous situations. We have seen 
in Bracken [4], Fricker [16] , Clemens 
[7] , Turkes[40], Lucas [29] that LSE 
method have been applied for eval- 
uating the parameters for fitting the 
homogeneous Lanchester equations to 
the historical battle data.The MLE 
method [33, 38] has not been explored 
particularly for fitting the historical 
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battle data till date. Also only one 
measure i.e. Sum-of-squared-residuals 
(SSR) has been explored for measur- 
ing the Goodness-of-fit (GOF). The 
main objective of this study is to as- 
sess the performance of the MLE ap- 
proach for fitting homogeneous as well 
as heterogeneous Lanchester equations 
to the Battle of Kursk. Various mea- 
sures of GOF [8] viz. Kolmogorov 
Smirnoff, Chi Square and R? have 
been computed for comparing the fits 
and to test how well the model fits 
the observed data. Applying the var- 
ious GOF measures considering the 
artillery strength and casualties of So- 
viet and German sides from the Kursk 
battle data of World War-II validates 


the performance of MLE technique.Section 
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2 presents in brief the 

overview of the battle of Kursk. Sec- 
tion 3 describes the mathematical for- 
mulation of likelihood estimation in 
case of both homogeneous as well as 
heterogeneous situations. Section 4de- 
scribes the Tank and Artillery data 
of Battle of Kursk and discusses the 
methodology for implementing the pro- 
posed as well as other approaches. 
Also, this section contains a perfor- 
mance appraisal of the MLE using 
various GOF measures. Section 5anal. 
yses the results after observing vari- 
ous tables and figures and discusses 
how well the MLE fits the data. Sec- 
tion 6 summarizes the important as- 
pects of the book.empirical 
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2 History of the 
Battle of Kursk 


After suffering a terrible defeat (With- 
drawal,) at Stalingrad in the winter 
of 1943, the Germans desperately wanted 
to regain the initiative. In the spring 
of 1943, the Eastern front was con- 
quered (Breakthrough) by a salient, 
200 km wide and 150 km deep, cen- 
tred on the city of Kursk. The Ger- 
mans planned in a classic pincer op- 
eration named Operation Citadel, to 
eliminate the salient and destroy the 
Soviet forces (Prepared Defense) in 
it. On 2 July 1943, Hitler declared, 
"This attack is of decisive importance 
and it must succeed, and it must do 
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so rapidly and convincingly. It must 
secure for us the initiative... The 
victory of Kursk must be a blazing 
torch to the world" [22, 24, 45]. 

The Germans started the Battle 
of Kursk on July 4, 1943 on the south- 
ern half of the Kursk salient, but this 
was merely to gain better artillery 
observation points(Hasty Defense). The 
battle began in earnest early in the 
morning of July 5, when Soviets con- 
ducted an artillery barrage before the 
Wehrmacht attacked. The Germans 
countered with their own planned bar- 
rages shortly thereafter and seized the 
initiative on both fronts. Soviet Gen- 
eral Rokossovsky redeployed his re- 
serves on the night of July 5 in order 
to attack the following day. The di- 
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visions of the 17th and 18th Guards 
Rifle Corps, with support from the 
3rd, 9th, 16th, and 19th Tank Corps, 
were beginning offensive operations 
at 5:30 a.m. on July 6 in support of 
13th Army. July 6 were considered 
as worst single day of CITADEL for 
German tank losses.On 7th July the 
German attack with armour forces in 
the northern and northeast side and 
captured the village Lutschki and con- 
tinued advancing towards the village 
of Tetrevino against the very strong 
Soviet infantry and armor battle. By 
evening July 7 the Germans were ableto 
capture the village Tetrevino. On 
8th July the German attacked with 
armour forces and barely captured 
the Teploye village. The Soviet coun- 
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terattacked and recapture the lost Teploye 
village.The soviet defended very well 
on that day although they lost 315 
tanks on that day in comparison to 
108 tank loss of the Germans. The 
German forces wanted to develop a 
sharp wedge towards Kursk via Oboyan 
village. The German forces attacked 
with more than 500 tanks.The soviet 
forces defended the Oboyan with so- 
phisticated artillery guns. Despite 
the strong defense the Germans were 
able to foothold over the Pena River. 
During the period of July 7-9both 
the sides had suffered largest num- 
ber of tank losses. Due to this reason 
the Germans planned to attack from 
less resistance Prokorovka side. Af- 
ter changing the direction of the at- 
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tack the German reached and seized 
the village of Novoselovka. The So- 
viets understood the German’s plan 
and they started using their reserve 
units. But despite of that the Ger- 
mans were able to break through (Break- 
through) the Soviet defenses by evening 
of the July 10. The German inten- 
tion at this point was to cross the 
river Psel to the extent of as many as 
troops and vehicles are possible. The 
Germans were able to seize a bridge. 
The engagement was at this point 
between artillery and tank.Both the 
forces were preparing for the battle 

of Prokhorovka. The Battle of Prokhorovka 
was the decisive phase of the Battle 

of Kursk. The Soviets started with 
artillery defense and later it turns 
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out to be totally tank against tank 
meeting engagement (Meeting Engage- 
ment,). The German tanks had to 


face minefield as well as well 


defended 


Soviet anti-tank weaponsDelaying, (De- 
laying,). The resulting titanic bat- 
tle was a tactical draw (Stalemate,). 
The Germans lost 98 tanks against 
414 Soviet tank losses. Hitler called 
off the battle. The losses from the 


fighting over July 12 and 13 


were ex- 


tensive on both sides. In KUTUZOV 


there was a heavy armour 


engage- 


ment between the two forces. The 
German forces destroyed 117 Soviet 
tanks. The Soviets also damaged 57 
German tanks.The battle on Proko- 
rovka still continued on 14th July. 


The German planned an offe 
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nsive op- 


eration named operation Roland. It 
started on July 14. The aim was to 
destroy the Soviet armor reservoirs 
(Deliberate Defense). The German 
armor units fought with the artillery 
forces that were defending the armor 
reservoirs of the Soviet in the south- 
ern part of the Prokorovka. Several 
tactical positions 

(Fortified Defense) were captured by 
the German forces. On this day the 
Germans were capable of performing 
minor offensive operations and they 
were launching attacks to form the 
Gostishchevo-Liski pocket. Hitler redi- 
rected it back to Isyum on July 16. 
According to the various war ana- 
lysts it is being considered that be- 
cause of Hitler’s decision, von Manstein 
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lost the availability of a powerful mo- 

bile formation that could have been 
very useful in the battle. During this 
time most of the engagements were 
between German infantry and Soviet 
tanks. Most of the damage was suf- 
fered by the Soviet because they were 

not equipped with modern antitank 
weapons that can deal effectively with 

the Soviet armour. The Soviets launched 
their counteroffensive along the Mius 
River on July 17. The Southwestern 
Front, commanded by Colonel Gen- 
eral Tolbukhin, attacked the heavily 
fortified Mius River line defenses(Fortified 
Defense). The Soviet counter attack 

was known as operation RUMANTSYEV. 
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3 Estimation 


Let S denote the time between two 
consecutive casualties for a side, its 
probability density function is denoted 
by fs(S). Let (mgt, ngt) represents 
the i” type force strengths (e.g. tank, 
artillery etc.) of blue and red forces 
of a battle for the k*” time instance 
respectively. Let us also denote (for 
k = 1,2... K) the time (a random 
variable) at which K'^ casualty oc- 
curs on Tj, (with realization tk). Let 
the Blue and Red casualties Xu and 
Y; in a combat are r.v. whose den- 
sities are defined by f., (s|a;, pi, qi) 
and fs, (s|b;, pi, qi) respectively where 
forms of the densities are known ex- 
cept the unknown parameters (aj, bi, pi, q;).1 
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is assumed that the of a random sam- 
ple (¿ga )from fs(S) can be ob- 
served. On the basis of the observed 
sample values (2; ,Y;, ) it is desired to 
estimate the value of the unknown 
parameters (a;,0;,p;,q;) . We fur- 
ther assume that the times between 
casualties are exponentially distributed, 
then the pdf of casualty for the Red 
(X) and Blue (Y) sides associated to 
the equation (1) and (2) can be rep- 
resented as in the equations 7 and 8: 
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—oo < di, pi, di < oo, Vi =1,2,...,F 


(7) 


F 


Ts (s) — TUARA 


i=l 


F 
ezp(—( ` bixi yu” )s) 
i=1 


—oo < bi, Pi, qi < oo, Vi A dl 
(8) 
The likelihood equation of n pairs of 
random variable (i;,9j;) is defined 
as the joint density of the n pairs 
of random variables, which is consid- 
ered to be a function of (a;i, bi, pi, qi). 
In particular, if (ży ,Y;,) are indepen- 
dently and identically distributed ran- 
dom sample from the density fs(S), 
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then the likelihood function is: . 
f(xy |ai, Pi, qi). f (yy lbi Di, Oi E 
f(xy |a; Di; qi)-F (u bi, Pi, qi) 
'Then, the joint pdf will be: 


L(a;,b i» Pis qi) TH (ai. Ti ¿Pia 4i) Èy 


i—li—l 


x (bia Pi du di y x 


ezp(— 63 aitu” yy 
biyi” xu t )8) 


To construct the likelihood func- 
tion from the available data set, it is 
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generally observed that casualty fig- 
ures are generally available at daily 
interval. Let L(a;,bi,pi,qi) be the 
likelihood function for the random vari- 
ables (¿y 95). If ai, bi, Di, ĝi are 
the values of a;, bi, pi, qi which maxi- 
mizes L(ai, bi, pi, qi) ñ then Qi, bi, Di, 
ĝi are the maximum-likelihood esti- 
mates of a;, bi, Pi, qi. Now, instead of 
maximizing the likelihood function we 
will maximize its logarithmic form 
since both the maximum values oc- 
cur at the same point and logarith- 
mic form is easily imputable. Thus, 
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on taking log of equation (8), we have 


F N 
lnL = UU] [eu (azia 2 )+ 


i=1t=1 


Yal (bixu yz P) — (aimi P yu t 


biyi” xu )s) 


Differentiating the Log-likelihood func- 
tion (9) partially with respect to a; 
and b; and equating it to zero, we 
have: 


and 


dinL Ya PON nod 
db; = b; — X ay?! yy ts =0 
t=1 i=1 t=1 (10) 
This gives 
N dy F N 
and 
N di F N 
> >. = Y on n (12) 
t=1 ^ 


Thus, the maximum likelihood esti- 
mates are: 


No. 
a; = >= Ti 


13) 
š F SN = =. Í 
Deiat tar Lip Pt Yu is 


43 


i 


No. 
2st Vi, 


F N : 
D Son yi, P Lyts 
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(14) 


4 Overview of the 
Kursk Database 


The Kursk Data Base (KDB) is de- 
veloped by Dupuy Institute (DPI) and 
is reformatted into a computerized 
database in1998. KDB is documented 
in the KOSAVE (Kursk Operation 
Simulation and Validation Exercise) 
[3]. The KDB contains daily on hand 
and losses for the four categories viz. 
manpower, tanks, APC and artillery 
for the Soviets and Germans for each 
of the 15 days of battle. Evidences 
of multiple force interaction in Kursk 
Battle shows multiple forces were fight- 
ing in the war. Therefore, developing 
heterogeneous model on this data is 
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justified. In the present study, we 
have considered only the tank and 
artillery data for developing hetero- 
geneous Lanchester model. Table 1 
shows the tank and artillery weapons 
on hand and losses during the 14 days 
of battle. Figure 1 shows a com- 
parison between the Soviet and Ger- 
man’s tank losses during the 14 days 
of battle. 


This book fits the generalized form 
of Heterogeneous Lanchester equations 
to the Battle of Kursk data using 
the method of Maximum Likelihood 
estimation and compares the perfor- 
mance of MLE with the techniques 
studied earlier such as the Sum of 
squared residuals (SSR), Linear re- 
gression and Newton-Raphson itera- 
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Figure 1: Comparison of daily num- 
ber of tank losses of the Battle of 
Kursk of WW II. 


tion.Different authors applied differ- 
ent methodologies for fitting Lanch- 
ester equations to the different battle 
data. The methodologies of Bracken, 
Fricker, and Clemen are applied to 
the Tank data of Battle of Kursk and 
results are shown in Table 3. 
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5 Goodness-of-fit Statis- 
tics 


First, we applied the technique of Least 
Square for estimating the parameters 

of the heterogeneous Lanchester model.The 
GRG algorithm [15, 31] is applied for 
maximizing the MLE and for min- 
imizing the LSE. For implementing 

the Least Square approach, the Sum 

of Squared Residuals (SSR) is mini- 
mized. The expression of SSR for the 
equation (1) and (2) is given as: 


14 


2 
SSR= >` (y — Y aim, P ys ty 
i=l 


t=1 
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14 2 
3 (g: = > biyy Pm). (15) 
t=1 i-1 

For implementing this expression from 
table 1 we have taken zero as ini- 
tial values for all the unknown pa- 
rameters. Then we start running the 
GRG algorithm iterative. The GRG 
algorithm is available with the Mi- 
crosoft Office Excel (2007) Solver [15] 
and MATLAB [31]. The GRG solver 
uses iterative numerical method. The 
derivatives (and Gradients) play a cru- 
cial role in GRG.We have run the 
program for 1000 iterations for get- 
ting the stabilized values of these pa- 
rameters. Once, we have the param- 
eters we compute the estimated casu- 
alties. With the difference between 
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the estimated and observed casual- 
ties, we computed the Sum of Squared 
Residuals. Similarly, we applied the 
GRG algorithm for optimizing the ob- 
jective function as given in equation 
(9).We check the graphs of estimated 
and observed casualties for both the 
LS and MLE based approaches and 
found that if we divide the data set 
into several subsets then we can im- 
prove the fit. As we increase the num- 
ber of divisions,the fit turns out to be 
better. The estimated casualty con- 
verges to the observed casualty. We 
have considered tank and artillery data 
for mixing the forces therefor eal (or 
b1) represents effectiveness of Soviet(or 
Germans) tanks against Germans(or 
Soviets) tanks and a2(or b2) repre- 
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sents effectiveness of Soviets (or Ger- 
mans)Arty against Germans (or So- 
viets) tanks.The variation of attri- 
tion rates throughout battle tells us 
how the different player in the battle 
performs. Whether they are acting 
defensively or offensively. 

The basic idea of using GRG al- 
gorithm is to quickly find optimal pa- 
rameters that maximize the log-likelihood. 
The objective is to find the parame- 
ters that maximize the log-likelihood 
or in other words provide the best fit. 
Given the values in Table 1, we in- 
vestigate what values of the param- 
eters best fit the data.Although we 
derived the estimates fora and b us- 
ing the MLE approach in equations 
(8) and (9), they are not applied di- 
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rectly. Log Likelihood is calculated 
using the equation (7) considering 0.5 
as the initial value of the parame- 
ters. Then, we optimized the entire 
duration of the battle of the likeli- 
hood function using the GRG algo- 
rithm.The model obtained after esti- 
mation of parameters is: 


tı =(1.46) ary 129, 404 
+ (.906)zy Syz 136 
ji =(.704) 1291 404 
+ (.953)y, 1985,13 


(16) 


As the data for the first day is 
extremely low, we drop it since it 
will pose a problem in the computa- 
tion of the likelihood and SSR func- 
tion. Also, the extremely low casualty 
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levels on the first day represent large 
outliers; thus, including the data of 
the first day affects the outcome to 
a great extent. Thus, the first day 
was dropped in fitting the data to 
the models. This approach is also 
justified by the historical account of 
the battle of Kursk, because the fight 
did not begin until July 5, the second 
day of the battle. Thus, dropping the 
data for the first day and dividing 
the remaining 14 days data into five 
phases, the total number of optimal 
parameters with each day as single 
phase is 102. This is a much bet- 
ter fit than any of the homogeneous 
model because both the residual as 


well as the likelihood are optimized.Log- 


likelihood is calculated using equa- 
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tion (7) and is maximized separately 
for each of the five phases. Let t 
denote the days, then the division 
is made as (to = t3), (ta = te), (tz — 
tg), (to = t11), and (tia = t15) Fitting 
the model over multiple phases re- 
sults in a better overall fit because 
there are additional parameters to ex- 
plain the variation in casualties. The 
model has been improved from par- 
titioning the battle into 14 phases. 
Each day of the battle is treated as 
mini-battle. For the purpose of com- 
paring models, R? value is calculated 
along with the Sum of squared resid- 


55 


uals (SSR). R2 value is calculated as: 


R?=1- — 


15; ^ 15 po 2 

t=1 (iy = ty) RI (Yi = 94) 
. PR 15 . me 

t=1 CH x ty)? 3:25:13 (g: sa yi)? 
(17) 


A larger R? value indicates better 

fit. Also, Goodness-of-fit measures 
namely; Kolmogorov-Smirnov statis- 

tic [8] and Chi-square ( x?) [8] have 

been calculated for the accuracy as- 
sessment of the MLE to that of the 
conventional approaches. Kolmogorov- 
Smirnov statistic is a measure of Goodness- 
of-fit, that is, the statistic tells us 

how well the model fits the observed 
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data. The Kolmogorov-Smirnov (KS) 
statistic is based on the largest verti- 
cal difference between the theoretical 
and empirical (data) increasing dis- 
tribution function. 


x. peser 
KS —maza «e eao[F (et) = ^30 ^ 
t* 2 
3 7 Fee) 
(18) 


where F(e;) is the cumulative distri- 
bution function of the estimated er- 
ror between the observed losses and 
the estimated losses for both sides. 
Chi-Square (x?) is another measure 
of Goodness-of-fit. Chi-Square is given 
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14 3 $ 
ease Gain? 
t=1 ii 
» P (19) 


(Yer — Yar)? 
T ^ 


t=1 Vi; 


where ti: and Èi are the observed 
and expected casualties respectively. 
This book also explores other his- 
torical battles and estimated goodness- 
of-fits statistics to analyse the mili- 
tary strategies for modeling the Lethal 
Behavior [2] of autonomous forces. 
These models are given in table 6. 
The directional field plots or D-field 
plots are a visual representation of 
the the system of differential equa- 
tions. These are shown in figure 7. 


58 


These graph gives an insight about 
the data and shows the divergent and 
convergent properties of the histori- 
cal battles.Generally we see that for 
developing combat model of aggre- 
gated forces different elements of the 
forces are aggregated together in terms 
of Combat Potential or Lethal Be- 
havior [2] and then differential equa- 
tions are used for representing their 
decay over time frame. These are 
pseudo-aggregated model, for the pur- 
pose of modeling heterogeneous forces 
are blended together and total aggre- 
gated losses are estimated and then 
again dis-aggregated for allocating these 
losses in different individual forces. 
For doing so a significant amount of 
information are being lost. In the 
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present approach these information 
loss is being managed by ignoring the 
aggregation and then dis-aggregation. 
Here we directly estimates the rela- 
tive attrition coefficients of one ele- 
ment against another one. The lim- 
itations of the current approach is 
that it does not consider other in- 
fluencing factor of a combat. If we 
follow the Dupuy’s QJMA concept 
along with the Lethal Behavior of au- 
tonomous forces [2, 13] for quantifi- 
cation and to define the Combat Po- 
tential and Lethal Behavior (CPLB) 
of an autonomous force 7 with n; au- 
tonomous fighting elements as: 


(20) 
mo; -pol; .poi.-Watt/Def 
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where, 


ni: number of autonomous ele- 
ments in the it” force 

OLI;: Operational Lethality In- 
dex of the it” element 

m4: mobility factor 

vi: vulnerability of the force 

l;: leadership 

t;: training 

e;: ethical 

mo;:morale [44] 

pol;: political 

po;: posture 

Watt/Def: weather effect on at- 
tacker or defender. 

The OLI of an autonomous mo- 
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bile weapon is defined as 


OLI Movite =((fmr) + p) 
RFA,AmC 
OLINon—Mobile =(r*tiRARıms. 
9CmdmhiAm)/d 
(21) 
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where 
f = firepower 
m = mobility 
r = radius 
p = punishment 
R = Rapidity 
F = Fire control effect 
A, = Ammunition Supply 
Am = Aircraft Mount 
C = Ceiling 
r* = Rate of fire 
t = Number of targets 
i = incapacitating 
R = range 
A = Accuracy 
Rı = Reliability 
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Ms = self-propelled mobility 
g = guidance 
Cm = charges(multiple) 
bm = barrel(multiple) 
h, = Wheel Track 


d = dispersion 


let us assume an autonomous Ar- 
mour Troop with 3 Tanks of Type 1, 
[hence Troop OLI is =Troop OLI(T'ank,) —3 
3.((fimairi) + pi) Ri FA As, Am, C1) 
is fighting with another autonomous troop 

Troop OLIrank,=3-((famer2) + 
p2)R2F2 As, Amo, C» 

Combat Potentials(CP) of these 
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sides are: 
CP, =n1 OLI, .m;.01.11.t1. 
mo; poi, Watt 
CP, —n».O LlI».mo».vo.l.t2. 


mo3.p02, Waef 


(22) 


Differentiating the above equations 
over time t we get, 


dC P. 
n Í =n1.OLI; .M1].U1 lat. 


(23) 


—n».O LlI».mo».v».l».t2. 


dWaeg 
dt 


dC P3 
t 


T1102. p02, 
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dividing above equations we have 


dC Pı NY AW att 
are 24 
ce wae M 


where c = OLI; .mi.v1.l1.t1,mo01.po1 
. ñ OLI2.m2.v2.l2.t2,Mm02.p02 ? 
if defender is at advantageous posi- 


2 dC P, 
tions then JCP; < 1, 


ë ni dW ate 
‘ng dWae f 


<1 (25) 


If weather effect is same for both the 


attacker and defender then Ae =1 
def 


it implies that 


13 .U3 .l5.t9, TILO5 .DO» OLI; .n; 


my,.v,.1,.t1, mMo,.po, OLI,.no 
(26) 
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5.1 Theorem-1: 

Art of Maneuvering 
Consider a combat between autonomous 
FORCE, (niJOLI,,m1,v1,l1,t1,mo1,p1) 
and FORCE) (noO LT», ma, V2, la, t5, TO», 
where mij, vi, li, ti, moi pi€ R 


CP, = — a.C P and CP, = —b.CP, 
(if direct fire mode) 
CP, = — a.C P, and CP, = —b.C P» 
(if indirect fire mode) 
(27) 


where a,b € R and for t = 0, CP, = 
CPi, and CP, = C P>, at t = 0 are 
real function, then: 


(i) if direct and indirect fire both 
are in cycle then adding above 
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two equations we have 
CP, = . +CPa) (28) 
and 


b 
CP, = -;(CPi + CPx) (29) 


If effectiveness of force is dif- 
ferent in direct fire mode and 
indirect fire mode then 


CP, = —a,.CP, = aə.C P> 
(30) 
and 


CP, = —b,.CP, — bə.C P> 
(31) 


(ii) Numerical strength increases 
to attacker whereas numerical 


68 


strength reduces to defender thus 
if side 1 is attacker and size 2 
is defender a tactical parame- 
ter is being introduced i.e. 


(33) 
CP, -adCP,  aCP,, 
Se (a). 
CP, pen. b CP, 


Therefore force ratio gets squared 
time increased. The ethical rea- 
soning of the autonomous forces 
are divided into two processes 
as we have seen in the Ethical 
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Governor in MissionLab[2]. These 
are Evidential Reasoning and 
Constraint Application. For con- 
strain application we follow the 
GOF mathematics as described 
above and for Evidential Rea- 
soning we adopted the Sun Tzu’s 
Art of War. Now let us con- 
sider Side 1 is attacker there- 
fore, 


CP, =n1.(d).OL1,.m,.01.11.t1. 
mo1.po1. Watt 

CP =n3.(1/d).OLI3.m3.02.l2.t2. 
mo» po» .W def 


(34) 


If the defender does not have 
the exact location of battle he 
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will spread across a large area 
which will reduce the force con- 
centration. Where as attacker 
is going to hit a particular point 
with much higher force concen- 
tration, therefore combat den- 
sity at hit point is: 


£ l4 
if CP, .dt.dl¡ and 
0 JO 
n (35) 
/ CP .dt.dlg 
0 J0 
where l5 > h 


Let us consider the defender is 
deployed and divided over 4 sub- 
units in 4 different positions, 
front, rear, right and left.Hence, 
CP» = OPE ront A Cheer 3 
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ORS OB = Ut OPE 
When time(t) and location(l) 
are unknown 


t 1 
| | CP,.dt.dl = 0 
0 0 


t l 
> n J Ut oC P5 dt.di = 0, 
0 ZO 
(36) 


Succor or Support 


P(t € T,l € L,Uj.g 


ae (37) 
J n C P) dt.dl < 0) 
0 0 


So even ng > n; but the suc- 
cour or support during an at- 
tack denoted by P(t € T,l € 


L, Uto fa Sa C P3.dt.dl < 0) so 
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the plan is (time,location) (t;, p;) 
is the governing factor, the like- 
lihood of the plan is 


L(t; €T, xi yi € X,Y) = 
le X,Y 
[i i CP! dzj.dy; < 0) 
(38) 


Soldier exceeds in number that 
shall matter nothing in the mat- 
ter of victory, consider a situa- 
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tion if n4 > ng and 


81 =P(t e T,l e L, Uo 


t 1 
J / C P! .dt.dn) 
0 0 


S; =P(t € T,l € L, Ufo 


t l 
f l C P$ .dt.dn2) 
0 JO 


(39) 


S} < 5? then S} will not con- 
tribute in victory SZ will con- 
tribute in victory. 


If enemy is stronger in number 
we can prevent him from fight- 
ing scheme is to know the plans 
and likelihood of success. The 
of a plan depends on the Com- 
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bat Potential (at location 1 at 
time t. therefore the likelihood 
of success of a plan is the ra- 
tio of combat densities of two 
sides. Therefore, 


L(x € win|ni, na, t, 0) = 
J: f c P.atal (40) 
J [' CPadtal 


Principle of activity i.e. the func- 
tional properties of CP as a func- 


tion of (lt) 

CP, (Lt) = ni (1,t).(d).O LI. 
m.4(1,t).u1(1, t).4.t1, 
mo poi (l, t). Ware (I, t) 

CP 2(1,t) = nə(1,t).(1/d).O LI. 
mə(1,t).uə (1, t).l2.t2, 
mo».po»(l, t) .Wae (1, t) 

(41) 

Taking logarithms on both sides 

we have 

logC P, (l, t) = log(n, (1, t)) + log(d) 

+ log(OLI,) + log(ma (l, t)) 

+ log vs (L2) + log(h) 

+ log(t1) + log(moi) 

+ log(por(1,t)) + log(Waw (1,8) 
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logC'Pa(1,t) = log(na(1,t)) + log(1/d) 

+ log(OLI,) + log(ma(l,t)) 

+ log(va(1,t)) + log(21) 

+ log(t2) + log(moz) 
( 


+ log(pox(l, t)) + log(Wae (l, t)) 
(43) 


Collect samples from various lo- 
cations C P(z!, y!,t), CP(z1 y!,t),.. 
C P(x",y".1") € Fy Carefully 
compare the opposing Army and 
Conceal Tactical Dispositions: 

that means we have to camou- 

flage the center location of the 

CP's and the distribution pat- 

tern or distribution dispersion. 
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5.2 Theorem-2: 
Art of Self Possession 


Victory depends not only on the 
n,m,l,mo,p,... other factors which 

is important is the force concentra- 
tion of CP at time t and location l. 
If comparatively it is more than the 
enemy then only victory can be pro- 
duced, it is not simply the number 
Ni, Vi; OLI;, li, ti, MOi, Pi, Watt/def it is 
combined effect over location and times. 


5.3 Theorem-3 
Art of Husbanding One’s 
strength 


The AOW of modern warfare is mainly 
evolved on three factors of conflict 
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dynamics, these are Goal in a strate- 
gic situation, time to reach the goal, 
potential to achieve the goal. In 
mathematical notation we define it 
as: 


1. Determine the goals (G1, Ga, ..., Gn) 


which are connected as network 
of systems like a Mosaic War- 
fare. 


2. Determine the time to reach near 
the Goal or at time when one 
wants to achieve the goal i.e. 
time (ta, ta, m tn). 


3. What is the potential to achieve 


the Goal G; at time t; i.e. CP(Gi,t;). 


The potential value is maximum 
at Goal G for simplicity we repre- 
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sent Goal as location l. The poten- 
tial value is gradually decreases ex- 
ponentially as goes far from the Goal. 
Consider system 


CPi, = n3.(d).OLI .m.vi.l. 
t1, mo, poi War 
C Pp = mə.(1/d).O LI,.ro.uo.lə. 


13, M02.p02.Wae f 
(44) 


Let Wa: be the function of l and t 
and War = sin(tx) with a,8 4 0, 
the CP value is maximum at loca- 
tion l and decays exponentially with 
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parameter estimated CP 


1 
B(Lt) = v/(Q.7.01.0| >) 


ezp(— 5 (Gn — DPE- ny) 


= CP 
(45) 
if oí = 02 = 0 
then 
1 
(2.7).0.] >. 1) 
= 1 (46) 


We can visualize this function as tac- 
tics and overall area under the curve 
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as strategies. so the CP(G;,t;) is not 
fixed and it varies with time. The 
overall distribution pattern of the tac- 
tics and strategies are represented by 
the theoretical distribution pattern 
of the equation (45).As we see in the 
curve all the situation or tactics are 
different from each other. But it can 
be imagined that situation 3 as well 
as situation 4 is mirror (inverted) im- 
age of situation 5 and so on. What 
remain constant is the middle situa- 
tion 1. Although these situations are 
imagined as ideal theoretical dimen- 
sion but in real practical time there 
will be not exactly with the theoret- 
ical dimensions . So in actual real 
condition the situation will be gov- 
erned by the Goal and time. So in a 
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practical sense we can analyse a ter- 
rain from the GIS and we can tag dif- 
ferent situation with time variation. 
That’s why San Tzu is imagining how 
the water flows on the ground . As 
water changes its direction and shape 
according to the terrain the potential 
army will have no fixed dimensions, 
it will change according to the ter- 
rain. So the efficient Commander has 
to decide three things to maintain his 
force as free flowing. Similarly this 
flow of water like force is dependent 
on the season also. Therefore in a 
particular time what should be the 
goal of today may not be important 
for tomorrow. Similarly the poten- 
tial to achieve the goal at time t; may 
not be equal to the required potential 
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to achieve the Goal at time t;. The 
art of husbanding One’s strength is 
first determine the correct Goal, and 
if the Goal is same to both the sides 
then one may try to achieve the goal 
before the opponents. And to achieve 
the goal before the opponent the po- 
tentiality also has to be increased or 
sufficiently large to well feed the force. 
This is the art of husbanding one’s 
strength. 

The density function is the tac- 
tics and its total area is the overall 
strategies. Do not repeat the tac- 
tics these tactics are governed by the 
plan is dependent on l and t. So 
a particular plan is just the realiza- 
tion during these period. So the i 
and £ should not be repeated 1 and 
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t are € L,R?. Military tactics are 
like Unto water. So strike on weak, 
Foe = Ground , Water = Soldier. 
As water has no constant shape war 
has no constant shape.Heaven born 
leader modify his tactics in relation 
to his desire. Earth, metal and planet 
are not important what is important 
is the seasonal changes of water a 
they one season gives way to the an- 
other season. this book uses the con- 
cept of Art of War [41] for establish- 
ing the above mathematical co nota- 
tions. 

The Matlab Code for generating 
this figure 
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eS |Loss ailes 5 
y=100*normpdf (x,0,1); 


3 pgen=polyshape(x,y); 


figure , plot (pgen) 


5 hold on 


xi-[-1.75:.1:1. 7515 

y1=100*normpdf (x1,0,1); 

pgen1=polyshape(x1,y1); 

plot (pgeni,’FaceColor’,’magenta 
2) 

hold on 

plot (polyshape([-1.75 -1.75 
1.75 1.75],[y1(1) 
,0.44,0.44,y1(1)]),? 
FaceColor”,*magenta?,? 
EdgeColor’,’magenta’) 

hold on 


E Se Coll 


y1=100*normpdf (x1,0,1); 


5 pgeni=polyshape(x1,y1); 
j| plot (pgeni ,’FaceColor’,’red’) 


hold on 
plot(polyshape([-.75 -.75 .75 
.75],[y1(1) ,0.44,0.44, y1(1) 
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19,*FaceColor","red',* 
EdgeColor’,’red’) 
hold on 
x= 952.1525] 
y1=100*normpdf (x1,0,1); 
pgeni-polyshape(x1,y1); 


3 plot (pgeni,’FaceColor’,’blue’) 


hold on 


s plot (polyshape([-.25 -.25 .25 


.25],[y1(1) ,0.44,0.44, y1 (1) 
10,*Facelolor”,*blue?, >? 
EdgeColor’,’blue’) 
hold on 
x1-[-1.25:.1:1.251;5 
yi-100*normpdf(x1,0,1); 
pgeni-polyshape(x1,y1); 
plot (pgeni,’FaceColor’,’cyan’) 
hold on 
plot (polyshape([-1.25 -1.25 
1.25 1.25],[y1(1) 
,0.44,0.44,y1(1)]),” 
Hace Color n cyan, 
EdgeColor’,’cyan’) 


33 hold on 
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40 


text(2,1,"Fortified Defence ( 
Hemmed-in)’,’Rotation’ ,90) 

text(1.5,1,’Deliberate Defence 
(Serious Ground)’,’Rotation 
290» 

text(1,1,’Delay (Open Ground)’, 
?^Rotation? ,90) 

text(.5,1,’Breakthrough (Facile 
Ground)’,’Rotation’ ,90) 

text(0,1,’Withdrawal ( 
Dispersive Ground)’,’ 
Rotation’ ,90) 

text(-.5,1,’Stalemate ( 
Contentious Ground)’,’ 
Rotation’ ,90) 

text(-1,1,’Hasty Defence ( 
Intersecting Highways)’,’ 
Rotation’ ,90) 

text (-1.5,1,’Prepared Defence ( 
Difficult Ground)”, * 
Rotation’ ,90) 

text(-2,1,’Meeting Engagement ( 
Desperate Ground)’,’ 
Rotation’ ,90) 
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13 label C Situations?) 
14 ylabel(’Strategy’) 
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5.4 Theorem-4 
Art of Studying Cir- 
cumstances 

Consider the battle 

CP, = af(CP,) + bg(CP;), 
CP; = af (CP2) + bg(CP;) 


(47) 


being a,b,c,d € R and f and g are 
smooth real function such that f(0) = 
g(0) = 0 then: 

(i) if abcd < 0 battle system (23) 
has no limit cycle. 

(ii) Assume that f and g are an- 
alytically 


f(CP,) = (CP)? + O(CP,)"and 
CP, = (CP,)**-! + O(C P)” 
(48) 
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for some positive integer number k 
and | where k Z l. Then there ex- 
its abcd such that system (23) has at 
least one limit cycle surrounded the 
origin which whenever exists is hy- 
perbolic. 

(iii) there exists f and g such that 
for same values of abcd system (23) 
has more than one limit cycle sur- 
rounding the origin. More over same 
values using g(C P2) = CP, that is 
for system (23). 
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5.5 Theorem-5 
Art of Studying Mood 


Consider the battle 
CP, = Af(CP,) + Bg(C P2), 


l (49) 
CP, = Cf (C P2) + Dg(CP1) 


being ABCD € R and iid random 
variables with N (0,1) Gaussian Dis- 
tribution and where f and g are smooth 
real function such that f(0) = g(0) = 
0 then the probability that it does 
not have periodic orbits is greater than 
or equalto 3. Equivalently, the prob- 
ability of having some limit cycles is 
smaller than or equal to i. then: 

(i) if abcd < 0 battle system (23) 
has no limit cycle. 

(ii) Assume that f and g are an- 
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alytically 
f(CP,) = (CPL)?! + O(CP,)*4and 


CP, = (CP) + O(C P2)” 
(50) 


for some positive integer number k 
and | where k Z l. Then there ex- 
its abcd such that system (23) has at 
least one limit cycle surrounded the 
origin which whenever exists is hy- 
perbolic. 

(iii) there exists f and g such that 
for same values of abcd system (23) 
has more than one limit cycle sur- 
rounding the origin. More over same 
values using g(C P2) = CP, that is 
for system (23) . 

the system does not have peri- 
odic orbits P(ABCD < 0) ABCD 
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and — ABCD have same distribution 
A,—A, P(ABCD < 0) = P(— ABCD < 
0) = P(ABCD > 0) Since P(ABCD = 
0) = OP(ABCD > 0) = P(-ABCD < 
0)= >, thus the probability of the sys- 
tem pairing atleast one limit cycle is 

< 3 


5.6 Theorem-6 
Art of Handling Large 
Masses 


Consider the random combat system 


CP, = A.f(CP,) + B.CP», 


(51) 
CP, = C.f (CPi) + D.CP,. 


where f(CP,) = a.CP,* + 
P fi( CP Y + B(CP,)™, with 
aß #0, k < 0 odd integers, m > 1 
and ABCD iid N(0,1) random vari- 
ables, Assume also that x=0 is the 
unique real root of f(C P.) 
= 0.Then: 

(i) when k > 1, the probability 
of having an odd number of limit cy- 
cles is 1/8, and the probability of not 
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having limit cycles or having an even 
number of there is 7/8. 

(ii) when k = 1 and 6 > 0, the 
probability of having an odd number 
of limit cycles is PF(a) < 1/2 and 
the probability of not having limit 
cycles or the probability of having an 
even number of them is 1 — P+ (a). 
Here P* : R — (0, 1/2) is a decreas- 
ing further that satisfies limy_,,, PT (a) = 
1/2, P+(0) = 1/8, limas% Pt (a) = 
0 an eae "pom 
Bre = a2 AIC )exp-^ E NN 
ae 
T(a) = (a,b, c, d) : ad — bc > 0; 
a(ao + d) € s when k = 1 and 8 < 
0, then some results as in item (ii) 
hold but changing P* by P^ where 
P*(o) = P*(—o) In all the cases , 
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each limit cycle is counted with its 
multiplicity. 
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6 Discussion 


Figures 2 and 3 show the graphs of 
Soviet and German Tank losses along 
with the losses estimated through max- 
imum likelihood approach. In this 
model a single set of parameters are 
estimated for representing the entire 
14 days of the battle.Figures 4 and 
5 show the performance of the same 
model when entire data set is divided 
into 5 phases. From these figures, it 
is apparent that fitting the models 
with division into 5 phases resulted 
in a much better fit.Figures 6 and 
7 show the further improvement in 
the data set by dividing it into 14 
phases where each day is considered 
as a mini battle. Further, the total 
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losses are divided into two compo- 
nents: Losses due to tank and Losses 
due to Artillery. The overall SSR 
and likelihood values are functions of 
pis and q;'s. Figures 8 and 9 shows 
the 3D surfaces and contour plots of 
SSR as a function of pı, q, and po, 
q» respectively.From these figures, we 
can see that the minimum SSR zone 
is represented by contours of 1.5 E--5 
and 2.5E 4- 5. Using a grid search 
in this zone, the best or optimal fit 
is obtained at pı = .129, qı = .404, 
p2 = .138, q2 = .136 with SSR 1.19E4- 
5. The ai, bi, as, ba values corre- 
sponding to the optimal fit are 

1.14, 0.70, 0.90, 0.95 respectively.Figures 
10 and 11 shows the surface and con- 
tour plots of likelihood as a function 
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of pi, qı and po, q2 respectively. From 
these figures, we can see that the zone 
of maximum likelihood is represented 
by contours of 6.0E + 3 and 5.0E + 
3 with MLE 5.11£ + 3.Using a grid 
search in this zone, the best or opti- 
mal fit is obtained at p, = .21, q1 = 
.28, pa = .02, qo — .04. The 01, bi, 
aə, bə values corresponding to the op- 
timal fit are 0.99, 0.88, 0.89, 0.96 re- 
spectively. Table 2 shows the results 
of Bracken [4], Fricker [16] , Clemen 
[7] and MLE approaches applied on 
the tank versus tank and artillery data 
under heterogeneous situation. This 
table shows the KS statistic for MLE 
(with 14 divisions) is 0.08674, which 
is less than any other estimation meth- 
ods implying that the method of MLE 


101 


fits better as compared to the other 
methods. Also, R2 is a measure of 
goodness of fit. Larger values of R? 
implies a good fit to the data. The 
R? value of MLE (with 14 divisions) 
isl. For comparing the efficiency of 
the different approaches, the root mean 
square error (RMSE) criteria is used. 
The RMSE of MLE with 5 divisions 
is 88.13 and the RMSE of MLE with 
14 divisions is .0005, which is found 
to be the minimum.The RMSE of Clemen’s 
Newton-Raphson Iteration model is 
116.19, which is found to be the max- 
imum. Therefore, efficiency(E) is mea- 
sured with respect to the RMSE of 
the MLE with 14 divisions. Thus, 
the E for MLE is maximum i.e. equal 
to 1 and E of Clemen’s model is min- 
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imum i.e. equal to 4.30E — 06. If the 
comparison is made among Bracken's, 
Fricker's and Clemens approaches, we 
can say that the Bracken approach is 
better. However, in all the cases the 
MLE outperforms other approaches. 
Based on all the GOF measures, it 
can be concluded that MLE provides 
better fits. In the present research 
we just demonstrate that if it is pos- 
sible for mixing two forces it is also 
applicable for more than two forces. 
'The number of parameters to be es- 
timated increases fourth folded for 
mixing one additional force. With 
the estimated parameters, we com- 
puted the casualty due to tank com- 
ponent and casualty due to artillery 
component (See Table 3).When the 
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14 days Battle data is considered with- 
out any division, a and b parame- 
ters are significantly small and a: > 
bı which implies German tanks were 
more effective than Soviet tanks.Similarly, 
when we compare a» against ba, ba > 
ag which implies Soviet artillery were 
more effective than German artillery. 
Table 4 shows maximised log-likelihood 
values with divisions into 14 phases 
where each day is treated as a mini- 
battle.Table 3 shows the optimal pa- 
rameters of heterogeneous Lanchester 
model with an R? of 1, RMSE of 
0.0005, chi-square of 1.9E-5, SSR of 
3.3E-6 and MLE of 13202. The pa- 
rameters are obtained from maximum 
likelihood estimation of heterogeneous 
Lanchester model of tank and Ar- 


104 


tillery data (table 1) from Kursk Database 
with each day as single phase. The 
GRG algorithm is applied for maxi- 
mizing the likelihood function given 

in equation (9).Also, the parameter 
estimates a;, bi, Pi, qi are given corre- 
sponding to the maximised log-likelihood 
values with divisions.From this table 

we can see that the patterns of the 
parameters for each day of the battle 

are same for both the sides. In addi- 
tion the tank component parameters 

are seen to be playing major role in 
the entire duration of the battle. Out 

of 14 days, 10 days the tank compo- 
nent parameters came out to be the 
maximum. That's why the result jus- 
tified the Battle of Kursk and was 
correctly termed as the largest tank 
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battle in the history. 
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7 Conclusions 


Although mathematical formulations 
are well established for heterogeneous 
Lanchester model, very few studies 
have been done to mode] actual bat- 
tle scenario. We have developed het- 
erogeneous Lanchester model for Kursk 
Battle from World War II using tank 
and artillery data. Allthe previous 
studies on Kursk Battle were done 
to capture the homogeneous weapon 
system (Tank against Tank or Ar- 
tillery against Artillery). The work- 
ing principles of this model were only 
applicable for homogeneous situation. 
So extending those models in hetero- 
geneous situation both theoretically 
and practically were main focus of 
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this book. We have formulated the 
likelihood expression under heteroge- 
neous situation and applied to fit model 
under heterogeneous Lanchester model 
for Kursk database.We have estimated 
the MLE of the different parameters 
that are proved to be statistically more 
accurate. The unfamiliarity to deal 
with the heterogeneous situation by 
the previous approaches motivated us 
to venture the minute details of the 
Kursk Battle. The estimates are cross- 
validated to control the problem of 
the over fitting. Also, these estimates 
possess the optimal properties of con- 
sistency, sufficiency and efficiency. So 
compared to the previous work, the 
present book opens up the opportu- 
nity for exploring the complicated struc- 
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ture of Kursk Battle of World War 
ILFigure: 8,Directional Field or D- 
Field Plot: Considerable insight into 
the Combat system described through 
equations as referred in the table 6 
can be gained by examining the be- 
haviour of the differential equations 
through the use of plots known as Di- 
rectional field plots or D-field plots. 
The graph gives an idea whether the 
system is going to stable at a point 
or diverges without actually knowing 
the solutions. These graphs show the 
convergent and divergent properties 
of the battles. 
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Figure 3 


: Fitted Losses plotted ver- 


sus real losses for the (a) German's 
Tanks Losses (b) Soviet's Tanks 
Losses without any division of the 
Battle of Kursk of WW II. 
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Figure 5: Total losses are divided 
into two components, losses due to 
tank and losses due to arty. 
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(b) Soviet's Tanks Losses 


Figure 6: Contour plot of (a) Log- 
likelihood and (b) SSR values for 
the tank data of Soviet and German 
sides of the Battle of Kursk of WW 
II. 
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Figure 7: 3D plot of (a) 
likelihood (b) SSR 
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Acronym and Glossary 


AI Artificial Intelligence. iv, v, 2 
AIC Akike Info Criterion. 18 
APC Armour Piercing Carrier. 16 


Artificial Intelligence Artificial in- 
telligence (AT) is the intelligence 
of machines or software, as op- 
posed to the intelligence of hu- 
mans or other animals. It is a 
field of study in computer sci- 
ence that develops and studies 
intelligent machines. Such ma- 
chines may be called Als.. v 


Autonomous Forces Autonomous 
forces are military systems that 
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can operate without human in- 
tervention. They are often used 


in unmanned vehicles, 


drones, 


and other military equipment. 
The Department of Defense has 
hypothesized that future demand 
for uncrewed systems will strain 


the capacity of the defense in- 


dustrial base.. vii 


BAIC Bozdogan’s Akike Info Con- 


sistency. 18 


Breakthrough A breakthrough oc- 
curs when an offensive force has 
broken or penetrated an oppo- 
nent’s defensive line, and rapidly 
exploits the gap.. 29, 33 


CHI Chi Square. 27 
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combat modeling Combat Model- 
ing focuses on the challenges in 
human, social, cultural, and be- 
havioral modeling such as the 
core processes of move, shoot, 
look, and communicate within 
a synthetic environment and also 
equips readers with the knowl- 
edge to fully understand the re- 
lated concepts and limitations.. 
iv 


CP Combat Potential. 64 


CPLB Combat Potential and Lethal 
Behaviour. 60 


DDPG Deep Deterministic Policy Gra- 
dient. 24 
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deep learning Deep learning is a type 
of machine learning that uses 
artificial neural networks to learn 
from data. Artificial neural net- 
works are inspired by the hu- 
man brain, and they can be used 
to solve a wide variety of prob- 
lems, including image recogni- 
tion, natural language process- 
ing, and speech recognition.. vii 


Delaying a defensive military action 
in which advance of an enemy 
is delayed by fighting as long 
as possible without the defen- 
sive force becoming involved in 
decisive battle that would im- 
peril its withdrawal. 34 


Deliberate Defense a defense or- 
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ganized before contact is made 
with the enemy and while time 
for organization is available; usu- 
ally includes a fortified zone (with 
pillboxes) and communication 
systems. 35 


DRL Deep Reinforcement Learning. 
24 


Fortified Defense A fortification (also 
called a fort, fortress, or stronghold) 
is a military construction de- 
signed for the defense of terri- 
tories in warfare, and is used to 
establish rule in a region dur- 
ing peacetime. This is also known 
as Dispersive ground according 
to AOW[37], with defPrep=90 
percent . 35, 36 
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GOF Goodness-of-Fit. 27 


GRG Generalized Reduced Gradi- 
ent. 11 


Hasty Defense Hasty Defense - this 
is a defense you assume when 
you are in immediate contact 
with the enemy. . 30 


KDB Kursk Data base. 45 

KS Kolmogorov Smirnoff. 27 

LSE Least Square Estimation. 23 

Meeting Engagement In warfare, 
a meeting engagement, or en- 


counter battle, is a combat ac- 
tion that occurs when a moving 
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force, incompletely deployed for 
battle, engages an enemy at an 
unexpected time and place. 34 


Military Strategy Military Strat- 
egy is the planning and execu- 
tion of the contest between groups 
of armed adversaries. It involves 
using all the military, economic, 
political, and other resources of 
a country to achieve the objec- 
tives of war. Military strategy 
is often divided into four com- 
ponents: ends (objectives), ways 
(courses of action), means (re- 
sources), and risk. Military strat- 
egy has a long history, dating 
back to ancient times. Some of 
the most influential thinkers on 
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military strategy include Sun 
Tzu, Machiavelli, Clausewitz, Jo- 
mini, Liddell Hart, and others. 
Military strategy is influenced 
by the nature of warfare, the 
characteristics of the adversaries, 
the goals of the policy, and the 
constraints of the environment. 
Military strategy is not static, 
but dynamic and adaptive to 
changing circumstances. vii 


ML Machine Learning. 3 


MLE Maximum Likelihood Estima- 
tion. 23 


Mosaic Warfare The Mosaic War- 
fare strategy is based on a con- 
ceptual framework that integrates 
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all battle platforms to establish 

a complete picture of a quick 

and decisive victory. It empha- 
sises the use of a system-of-systems 
network to achieve this goal. This 
approach is applicable for multi- 
spectrum conflict(army, navy , 
airforce, cyber,information, AT) [28]. 
79 


OLI Operational Lethality Index. 61 


Prepared Defense preparedness re- 
quires considerable and consis- 
tent conceptual clarity, facile ground (q 
29 


QJMA Quantified Judgement Mod- 
eling Analysis. 60 
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RL Reinforcement Learning. 26 


rlFiniteSetSpace A discrete obser- 
vations or actions space in Re- 
inforcement Learning. 26 


rlNumericSpace A Continuous ob- 
servations or actions space in 
Reinforcement Learning. 26 


RMSE Root Mean Square Error. 102 
SSR Sum of Square Residual. 16, 
46 


Stalemate Stalemate is a situation 
in chess where the player whose 
turn it is to move is not in check 
and has no legal move.. 34 
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Sun Tzu’s Art of War The Art of 
War is an ancient Chinese mil- 
itary treatise that is attributed 
to Sun Tzu, a legendary mili- 
tary strategist and philosopher. 
The book consists of 13 chap- 
ters, each focusing on a differ- 
ent aspect of warfare, such as 
strategy, tactics, terrain, intel- 
ligence, and leadership. The 
book is widely regarded as one 
of the most influential works on 
military theory and practice in 
history, and has been applied 
to various fields beyond war- 
fare, such as business, politics, 
sports, and culture. vii 


Withdrawal A tactical withdrawal 
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or retreating defensive action is 
a type of military operation, gen- 
erally meaning that retreating 
forces draw back while main- 
taining contact with the enemy. 

. 29 
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Index 


adaptive, 7 

Aggregation, 8 

ambush, 15 

Ardennes Campaign, 
15 

attrition, 6 

autonomous forces, 
2 


Battle of Kursk, 4 
Bayes, 11 


Chi-square, 5 
combat modeling, iv 


Deep learning 
Neural 
Net- 
work, 
12 

differential equa- 
tion, iv 


estimation, 11 


generalized reduced 
gradi- 
ent, 11 

Goodness-of-fit, 5 


heuristic, 7 

High Resolution 
Simu- 
lation 
Model, 
8 


Killer Robots, 6 


Kolmogorov- 
Smirnov, 
5, 56, 
57 
Kursk, 3 
Lanchester’s 
Heteroge- 
neous 
Model, 
5 
Homogeneous 
Model, 
5 


linear law, 12 

Lanchester’s 
Equation, 7 
square law, 


12 
learning, 2 
reinforcement 
learn- 
ing, 3 


Least Square, 11 
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likelihood of success, 


74 


Maximum Likeli- 
hood, 
11 

Method of Moments, 
11 

Mosaic Warfare, 79 


Sum-of-square resid- 
uals, 5 
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supervised, 3 


theater model, 7 


unsupervised, 3 


weapons of mass de- 
struc- 
tion, 7 

withdraw, 14 


